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We quantify a limitation in the usual accounting of the finite-size effects, where the leading [(Za)'*] and 
subleading [{Za)^\ contributions to the Lamb shift are given by the mean-square radius and the third Zemach 
moment of the charge distribution. In the presence of any non-smooth behaviour of the nuclear form factor at 
scales comparable to the inverse Bohr radius, the expansion of the Lamb shift in the moments breaks down. This 
is relevant for some of the explanations of the “proton size puzzle”. We find, for instance, that the de Riijula 
toy model of the proton form factor does not resolve the puzzle as claimed, despite the large value of the third 
Zemach moment. Without relying on the radii expansion, we show how tiny, milli-percent (pcm) changes in the 
proton electric form factor at a MeV scale would be able to explain the puzzle. It shows that one needs to know 
all the soft contributions to proton electric form factor to pcm accuracy for a precision extraction of the proton 
charge radius from atomic Lamb shifts. 


I. INTRODUCTION 

The proton structure is long-known to affect the hydrogen 
spectrum, predominantly by an upward shift of the ^-levels 
expressed in terms of the root-mean-square (rms) radius. 

Re = s/W}e, {r^)E = Jdrr^PEif^, ( 1 ) 


A Lorentz-invariant definition of the above moments is 
given in terms of the electric form factor (FF), Ge{Q^), as: 


{r^)E = 
{r^)Ei2) = 



GeIQ^) 
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of the proton charge distribution pE- leading order (LO) in 
the fine-structure constant a, the nth S'-level is shifted by (cf., 
[ 1 ]): 


af;„s(lo) 
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( 2 ) 


where Z = 1 for the proton, irir is the reduced mass. The 
proton charge radius has thus been extracted from the hydro¬ 
gen (eH) and muonic-hydrogen (/rH) Lamb shifts, with rather 
contradictory results: 


REpieU) = 0.8758(77) fm [2], (3a) 

REpip-H) = 0.84087(39) fm [3,4]. (3b) 

The eH value is backed up by the extractions from electron- 
proton (ep) scattering [5, 6], albeit with a notable excep¬ 
tion [7]. 

The next-to-leading order (NLO) effect of the nuclear 
charge distribution is given by [8]: 


At the current level of precision, the eH Lamb shift sees 
only the LO term, while in pH the NLO term becomes appre¬ 
ciable. An immediate resolution of the eH vs. pH discrepancy 
(aka, the proton size puzzle) was suggested by de Rujula [10], 
whose toy model for proton charge distribution yielded a large 
Friar radius, capable of providing the observed pH Lamb shift 
using the Re value from eH. Shortly after, this model was 
shown to be incompatible with the empirical FF Gb extracted 
from ep scattering [11, 12]. In this work we find that the pH 
Lamb shift in de Riijula’s model is not described correctly by 
the standard formulae of Eqs. (2) and (4). The correct result 
involves an infinite series of moments, and it does not provide 
any significant reduction of the discrepancy in that model. We 
shall consider a different scenario of mending the discrepancy 
by a small change in the proton FF, using the corrected for¬ 
mulae. 


II. LAMB SHIFT: TO EXPAND OR NOT 


Ai?„s(NLO) = 

with Re( 2 ) = ^{f^)E{ 2 ) the Friar radius and 


( 4 ) 


{r^)Ei 2 ) = j drpE{r) J dr' \r - r'\^ PE{r') (5) 


the third Zemach moment. Other effects of proton struc¬ 
ture, such as polarizabilities, play a lesser role in both normal 
and muonic hydrogen, and are not in anyway of relevance to 
the present discussion of finite-size effects. 


Our main observation is that the standard expansion in the 
moments is only valid provided the convergence radius of the 
Taylor expansion of Ge in is much larger than the in¬ 
verse Bohr radius of the given hydrogen-like system. In other 
words, for ~ (Zarrir)'^, the electric FF must be repre¬ 
sentable by a quickly convergent power series. 

To see this we write the electric FF correction to the 
Coulomb potential (—Zajr) as follows: 

VMr) = — ImGEit), (7) 
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FIG. 1: One-photon exchange graph with FF dependent electromag¬ 
netic vertex. 


where Im Ge is the discontinuity in the FF across the branch 
cuts in the time-like region. This potential is derived by taking 
the non-relativistic limit of the one-photon exchange graph in 
Fig.l, and assuming the FF to satisfy the subtracted dispersion 
relation: 


with 

- 1 + . ( 12 ) 

does not exactly reproduce the Friar term in Eq. (4). To do so, 
we need to treat the potential to second order in perturbation 
theory. The second-order Lamb-shift is, to 0{a^), given by: 
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Here, as usual: = —q^ > 0, with q the four-momentum 

of the exchanged photon; to > 0 is the particle-production 
threshold associated with the given discontinuity of the FF. 
The use of the dispersion relation here is very much analo¬ 
gous to the Schwinger’s method of calculating the vacuum 
polarization (Uehling) effect [13]. 

The effect of this Yukawa-type correction on the hydrogen 
spectrum can easily be worked out using time-independent 
perturbation theory. At first order we obtain the following 
expression for the 2P-2S Lamb shift: 
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FF(1) 
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Combining it with the first order, we see that the third moment 
is replaced by the third Zemach moment, leading to Eq. (4). 
Certainly, if the expansion of the first-order result in moments 
breaks down, this second-order calculation is inadequate too. 
In either case, however, the second-order shift is suppressed 
by Za, with respect to the the first-order shift. 

It is furthermore convenient to express the exact (unex¬ 
panded) first-order Lamb shift in terms of Ge, rather than of 
its discontinuity. Observing that in any identity of the form, 

-| pOO pOO 

- dtW{t)lmGEit) = dQwiQ)GE{Q^), (14) 
Jto Jo 

the function W is the Stieltjes integral transform of w, i.e.: 


W{t) = 



w{Q) 


(15) 


we recast Eq. (9a) in terms of Ge by computing the inverse 
Stieltjes transform [14] of 


where in the second step we have expanded in the moments 
of the charge distribution, introduced in Eq. (1), using the fol¬ 
lowing (Lorentz-invariant) definition: 


{rnE = 


(N + iy. 



Im GE{t) 

fN/2+l • 


( 10 ) 


One can now see that the convergence radius of the power- 
series expansion in moments is limited by to, i.e. the proxim¬ 
ity of the nearest particle-production threshold. For the pro¬ 
ton we of course expect to to be a hadronic scale, of which 
the pion mass is the lowest, and the series should converge 
quickly for hydrogen, and in fact for most of the hydrogen¬ 
like systems. The new physics search is a different matter. 
Inclusion of new light particles, or anything that changes the 
FF at very low (and equivalently, the charge distribution at 
large distances) may invalidate the expansion in the moments. 
Before illustrating this point more quantitatively (in Sect. Ill), 
we need to make a few more technical steps. 

First, we observe that the expansion of the first-order Lamb 
shift to order a®. 


,FF(1) 
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We also can express the Lamb shift in terms of the 
spherically-symmetric charge distribution, knowing that the 
latter is given by the following Laplace transform of the FF 
discontinuity: 


1 r°° 

PeY) = / df ImG'£;(f)e""''^. (18) 

{2TrY r Jt„ 


Collecting it all together, we have 
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FIG. 2: Integrand of the Ist-order contribution to the Lamb shift [cf. 
Eq. (I9a) with Eq. (17)] of eH (blue dash-dotted line) and of /rH (red 
solid line) for the dipole FF. The dotted vertical lines indicate the 
inverse Bohr radius of the two hydrogens, while the vertical dashed 
line indicates the onset of the ep data. 


with w{Q) given in Eq. (17), and the Bohr radius a = 
l/{Zamr). The expression in terms of the charge density 
is the simplest and has the most intuitive interpretation: the 
hrst-order Lamb shift is given by the mean-square radius cut 
off at the Bohr radius by the Coulomb wave function. Indeed, 
it is simply the LO result, Eq. (2), with replaced by 

The expression in terms of the EE, Eq. (19a), can be derived 
without the reference to the dispersion relation. Instead, one 
can follow the momentum-space derivation of the Wichmann- 
Kroll contribution in Ref. [15]. The weighting function w{Q) 
then arises as a convolution of the hydrogen wave functions, 
which is interpreted as the atomic EE. The use of the disper¬ 
sion relation for the proton form factor will, however, allow 
for a simple and systematic account of the retardation effects, 
as will be discussed elsewhere. 

We are now in position to plug in a specific charge density 
or a EE parametrisation and see how good the expansion in 
the moments is. Eor instance, we hnd that the expansion does 
not converge for the charge density proposed by de Rujula to 
explain the proton size puzzle. His model does give an un¬ 
usually large Eriar radius, however it does not affect the /iH 
Lamb shift in the way expected from the expansion into mo¬ 
ments. We hnd in this model that, in the corner of the param¬ 
eter space used to explain the puzzle, the /rH Lamb shift ob¬ 
tained using Eq. (19) is very similar to the LO result, Eq. (2), 
and is nearly not as sensitive to the Friar-radius value as the 
NLO term, Eq. (4), by itself. In contrast, for most of the pop¬ 
ular FF parametrisations used in htting the ep data, including 
the dipole form, the expansion in moments is appropriate. We 
have checked all of the parametrisation listed in Appendix of 
Ref [16]. 

In general, the breakdown of the moment expansion can 
only be expected for FF with non-smooth behaviour in the 
region of Q comparable to inverse Bohr radius. To see this 
more explicitly we plot the integrand of Eq. (19a) in Fig. 2, 
for the dipole form of the proton FF, i.e.: Gep = (1 + 


(5^/0.71GeV^)“^. One can see that the total integral over 
Q is the result of large cancellations around the Bohr radius 
scale. Any relatively small variation in the FF around that 
scale may lead to significant effects in the Lamb shift. In 
order to define a proper charge radius which could be deter¬ 
mined from both the atomic and scattering experiments it is 
therefore mandatory to decompose the FF into the “smooth” 
and “non-smooth” parts. The contribution of the former can 
then be expanded in moments, while the latter must be treated 
exactly. This is already done in the case of QED contributions 
to the FF. 

Concerning the proton size puzzle, let us note that the two 
curves in Fig. 2 represent the two hydrogens, and that their 
regions of large cancellations are well separated. Also sep¬ 
arated is the region of the existing ep scattering data, > 
0.004GeV^. Obviously, if there is a small missing effect in 
the FF responsible for the puzzle, it must be localised near 
one of the two inverse Bohr radii, where its impact is max¬ 
imised. It also should not contradict the ep data, hence the 
allowed region is anywhere from 0 to roughly 50 MeV. On the 
other hand, only in pH the radius is extracted from the classic 
Lamb shift, the extraction from H combines a number of dif¬ 
ferent transitions. Hence, the simplest way to solve the puzzle 
is to have the missing effect in the pH Lamb shift, i.e., to place 
it between 1 and 50 MeV. To see this more quantitatively, we 
next consider a toy model for the possible effect in the FF in 
that region. 


III. RESOLVING THE PUZZLE 


We assume the electric FF to separate into a smooth {Ge) 
and non-smooth part (Ge), such that, 

Ge(Q2) = Gs(Q") + Gb(Q2). (20) 


For the smooth part we shall take a well-known parametrisa¬ 
tion which hts the ep data, while for the non-smooth one we 
consider a single Lorentzian:^ 


Ge(Q^) = - 

TT 





( 21 ) 


which describes a Breit-Wigner type of peak around Qq with 
the width given by 2e^, and which does not contribute to the 
charge: G_e(0) = 0. 

According to Fig. 2, in order to make a maximal impact on 
the puzzle, the fluctuation Ge must be located at the extremi 
[of w{Q) in Eq. (19a)] around of either eH or pH inverse Bohr 
radius. Here we shall only consider the latter case and set: 


Qo =+ o^lM657MeV. (22) 


^ This anzatz turned out to be similar in spirit to the one considered by Wu 
and Kao [17], who, however, deduced its effect on the Lamb shift through 
the expansion in moments, thus leading to incorrect results. 
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£ [MeV] 



Eq. 

Ge 

Ge 

Ge 

(r^)B [fm^] 

(r^}E [fm^j 

(6a) 

(12) 

(0.9014)2 

(1.052)^ 

-(0.2016)^ 

(6.384)® 

(0.8786)® 

(6.394)® 

Lamb-shift, expanded 
<P-2s(eH)[neV] 
752P-2s(FH)[peV] 

(11) 

-0.65690 

-4202 

0.03684 

4961 

-0.62006 

759 

Lamb-shift, exact 
^^2P-2‘?(eH)[neV] 
i?2pi\’s(FH)[peV] 

(19a) 

-0.65691 

-4202 

0.03684 

551 

-0.62007 

-3651 


TABLE I: Lamb shift and moments corresponding to our model FF, 
with e = 0.253 MeV, A = 5.65 x 10"® MeV^ and Qo set by 
Eq. (22). 


FF is extremely tiny. 


FIG. 3: Parameters of Gb for which the eH and Lamb shifts of 
Eq. (24) are reproduced, given Qo set by Eq. (22). 


This means we place a small fluctuation exactly on top of the 
maximum above l/a^n in Fig. 2. The sensitivity of our re¬ 
sults to the precise value of Qq is very mild for Qo between 1 
and 10 MeV. 

This choice of Qo > ^/a^H conditions the choice of the 
smooth part, in case one wants to solve the puzzle. Indeed, 
since the non-smooth part affects mostly the pH result, the 
smooth part must have a radius consistent with the eH value. 
We therefore adopt the chain-fraction fit of Arrington and Sick 
[18]: 


Ge{Q^) 


1 

3.478Q:^ ~ 

0.140 
1.311 Q'^ 

' ~ 1.128 

1 - 0.233 


(23) 


The other two parameters parametrising Ge, A and e, are 
fixed by requiring our FF to yield the empirical Lamb shift 
contribution, in both normal and muonic hydrogen, i.e.: 

■^ 2 .piTs^(eH) = -0.620(11) neV, (24a) 

= -3650(2) qeV. (24b) 

Note that these are not the experimental Lamb shifts, but 
only the finite-size contributions, described by Eqs. (2), (4), 
with the corresponding empirical values for the radii. In 
the eH case we have taken the CODATA value of the pro¬ 
ton radius, Eq. (3a), which is obtained as an weighted aver¬ 
age over several hydrogen spectroscopy measurements, and 
Re( 2 ) = 2.78(14) fm [9]. In the pH case we have taken the 
values from Ref. [4], hence Eq. (3b) for the radius and the 
same as the above value for Re{ 2 )- 

Eigure 3 shows at which A and e our EE complies with ei¬ 
ther eH (dot-dashed curve) or pH (solid curve) Lamb shift. 
Eor A = 5.65 x 10"® MeV^ and e = 0.253 MeV, our EE 
describes them both, thus solving the puzzle. 

We emphasise that the magnitude of the fluctuation in the 


\Ge/Ge\ <3x10-^, (25) 

for any Q^. Nevertheless, it obviously has a profound effect 
on pH Lamb shift. Its effect on the second and third moments 
is given by: 


(r^)E 

P)e 


12iA f 1 1 


(26) 


^ I (-Qo - (-Qg+ie2)5/2 
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+ 0(eVQo)- 


(27) 


The numerical values of these moments, together with their 
“would be” effect on the Lamb shift and the non-expanded 
Lamb result, are given in Table I. One can see that the expan¬ 
sion in moments breaks down for the fluctuating contribution 
to pH. 

We furthermore checked that by moving the bump (i.e., the 
Lorentzian structure) in to higher Q, anywhere between 
1 and 50 MeV, we can obtain a stable solution conforming 
both hydrogens. However, starting from Qo « 5 MeV the 
width and the magnitude of the bump quickly increase, reach¬ 
ing percent effects for Qo > 20 MeV. We find thus that the 
existing ep data, which start at about 60 MeV, and have the 
precision of a few percent, rule out the existence of such a 
“solving-the-puzzle bump” anywhere above 50 MeV. New ep 
experiments at lower Q of similar precision could lower this 
exclusion limit to about 20 MeV. 

To summarise, introducing a tiny (less than 30 ppm) bump 
on an empirical EE used to fit the ep scattering data, we are 
able to reconcile the eH and pH Lamb shift results. The bump 
is localised near the inverse Bohr radius of pH, and as such 
affects mostly the pH result. It has a small effect on the charge 
radius and the eH Lamb shift, and certainly no effect on the fit 
of the present ep data. 

Of course, until one finds a physical justification for such a 
structure at low Q, one should take this model only as a word 
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of caution against the very optimistic view of uncertainties 
in the charge radius extractions. All the contributions to the 
nucleon (as well as the lepton) electric FF at low Q should be 
known to better than a percent accuracy before claiming the 
discrepancy. 

IV. CONCLUSION 

We have shown that the effect of the nuclear charge distri¬ 
bution on the hydrogenic Lamb shift is not always expand¬ 
able in the moments of charge distribution — a very small 
fluctuation in the charge distribution near the Bohr radius of 
the atomic system may invalidate the expansion. The exact 
(non-expanded) expression for the finite-size effect is given 
by Eq. (9a), or equivalently Eq. (19). 

An example of charge density which invalidates the expan¬ 
sion of the pH Lamb shift in the moments is provided by the 
toy model of de Rujula [10]. Applying Eq. (19) to this model 
shows no significant reduction of the discrepancy between the 
eH and pH extraction of the proton charge radius. 


We have set up another simple model for the proton electric 
form factor which resolves the puzzle. Although we do not 
insist that this model has anything to do with reality, it demon¬ 
strates how a tiny (« 10“®) effect on the proton form factor, 
localised at low Q, may invalidate the usual line of arguments 
leading to the discrepancy between electronic and muonic ex¬ 
periments. We hope it will motivate a theoretical search of 
the possible low-Q effects in the form factors, as well as the 
new experiments called to constrain or exclude this scenario 
of resolving the proton size puzzle. 
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